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I. INTRODUCTION. 

A common practical requirement for a nonlinear sys¬ 
tem to have localized solutions is the stability of the zero 
background. This requirement on the one hand ensures 
absence of growing small fluctuations far from the non¬ 
linear mode. More importantly, this means the stable 
existence of a system itself without any excitations, i.e., 
in the “vacuum state”, in which the system is prepared 
experimentally and against which nonlinear modes are 
excited. In the conservative case this implies possibility 
for stable propagation of linear waves (real eigenvalues of 
the linear Hamiltonian), while in dissipative systems this 
means decay of all small amplitude excitations (the back¬ 
ground is an attractor with a nonzero basin). In this con¬ 
text parity (P) - time (T) symmetric [Ei systems repre¬ 
sent a special case which on the one hand obey gains and 
losses, and one the other hand may have pure real spectra 
in some domains of the parameter spaces [2| (this situa¬ 
tion is referred to as unbroken PT-symmetric phase Q) 
allowing for propagation of linear waves and possessing 
continuous families of solutions. Therefore linear poten¬ 
tials like parabolic Q, Scarff II Q, or PT-symmetric ex¬ 
tension of Rosen-Morse II potential [5[ obeying pure real 
spectra and supporting localized modes, received partic¬ 
ular attention also from the point of view of the existence 
of nonlinear families; see Ref. [6[, Refs. ( 3 , 111 , 13 ? 

respectively. 

Nonlinear modes can also be found in a region where 
the linear PT-symmetry is broken [ 3 ]. Moreover, fam¬ 
ilies of nonlinear modes can be stable ESini in those 
regions. However, the mathematical existence of such 
nonlinear modes, and of even stable ones, does not yet 
guarantee their practical utility, because the way of their 
excitation in a system which is linearly unstable in its 
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“vacuum” state, remains questionable. This leads us to 
the first goal of the present paper, which is a suggestion 
on how stable nonlinear modes ean he exeited in systems 
where the linear VT-symmetry is broken. The idea is 
based on the possibility of “switching-on” nonlinearity si¬ 
multaneously with gain and dissipation. Such possibility 
can be implemented, in particular, when the nonlinearity 
and gain-and-loss strength are characterized by a single 
parameter (let us call it e) and disappear when this pa¬ 
rameters becomes zero (e = 0). If at e = 0 the system is 
Hamiltonian, it allows for stable propagation of the lin¬ 
ear modes, and the only stability issue which has to be 
verified is the stability of the solution branch e > 0, bi¬ 
furcating from e = 0. Then, if the stability is confirmed, 
one can consider the adiabatic grows of e = e(t) in time 
as a way to excite a nonlinear mode, which persists stable 
even if the final value of e corresponds to the underline 
linear system with broken PT-symmetric phase. 

To solve this problem mathematically, a suitable 
framework is the use of potentials supporting some exact 
solutions. Such complex potentials can be constructed, 
say using the “inverse engineering” as this was suggested 
in (here we also mention other examples of particular 
exact solutions for PT-symmetric potentials published 
in Bill). The method consists in computing a complex 
potential starting with an a priori given solution, which 
however must satisfy some constraints to ensure the exis¬ 
tence of the potential. Generally speaking, the potentials 
supporting particular exact solutions appear of rather so¬ 
phisticated forms, what may constitute a significant dif¬ 
ficulty for their practical implementation. Therefore as 
a complementary task of this paper we consider genera¬ 
tion of the modes in potentials of relatively simple and 
experimentally feasible forms, bearing in mind their ap¬ 
plications in optics of atomic gases [IJ well as in the 
PT-symmetric physics of Bose-Einstein condensates [Hl- 
113, where PT-symmetric potentials can be created and 
modified in situ. 

The rest of the paper is organized as follows. In Sec.HIl 
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we describe the general theory and approach to the non¬ 
linear Schrodinger equation with the PT-symmetric po¬ 
tentials allowing for exact particular solutions. In Sec. Hill 
we discuss applications of the method to the parabolic 
potential and in Secs. IIVI and m we consider linear and 
nonlinear modes in both single-well and double-well po¬ 
tentials, respectively. Particularly, the problem of the 
nonlinear mode excitations in the single-well and double¬ 
well T^T-symmetric potentials allowing for particular ex¬ 
act solutions is also addressed in Secs. HYi and El Our 
results are summarized in the Conclusion. 


II. GENERAL THEORY AND APPROACH 

We consider the nonlinear Schrodinger (NLS) equation 
with the T^T-symmetric potential and space-modulated 
nonlinearity (abbreviated below as T^T-NLSE) 

idtip = -^dlip + Ue{x)'ip + Ge{x)\i>\'^'ip, ( 1 ) 

where dt = d/dt^ dx = d/dx^ ^ is the com¬ 

plex envelope of the electrical field, Ue{x) = U*{—x) (i.e., 
Re[U,{x)] = Re[Ue{-x)], lm[Ue{x)] = -lm[U,{-x)]) and 
Ge{x) describe the complex-valued linear T^T-symmetric 
potential and real-valued inhomogeneous nonlinearity, re¬ 
spectively, and an asterisk stands for the complex con¬ 
jugation. For specific the real-valued potential, i.e., if 
U^{x) = U*{x)^ Eq. ^ reduces to the conservative NLS 
equation with space-modulated linear and nonlineari- 
ties coefficients. In that case, particular exact solutions 
and corresponding dynamical behaviors were extensively 
studied in literatures (see, e.g., [2ll-[^b 

In order to implement the procedure described in the 
Introduction we assume that: 
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U,{x) = Y, + ieW{x), (2) 

j=o 

Ge{x) = e'^Gix), (3) 

where Vj{x){j = 0,1,2) are constituents of the real 
potential, W{x) is the real gain-and-loss distribution, 
e > 0 is the bifurcation parameter parameterizing a 
branch of the solutions: we emphasize that generally 
speaking e is not considered small. To ensure the VT- 
symmetry we consider Vj{x) = Vj{—x) {j = 0,1,2) and 
W{x) = —W{—x). At e = 0, Eq. ([T]) becomes the linear 
Schrodinger equation. Our main interest will be focused 
on the case where W{x) ^0 (i.e. when the linear poten¬ 
tial Ue{x) is non-Hermitian). 

We concentrate on stationary solutions of T^T-NLSE 
m in the form '0(x,t) = (f){x)e where /r is a real 
spectral parameter, and the complex valued nonlinear 
eigenmode (j){x) satisfies the stationary NLS equation 

= + Ue{x)4) + Ge{x)\4)\'^4>, (4) 


subject to the zero boundary conditions lim^^^ioo = 

0 . 

Now, following the strategy described in the Introduc¬ 
tion we assume that the solution of the linear eigenvalue 
problem 

I (p 

-^O0n(^) ~ -) L/q = “ ^^2 (5) 

where n = 0,1,2,..., is known. In Eq. ([5]) we assume 
also that the spectrum is discrete (what is applied for 
all examples considered in this paper) and distinguished 
the eigenfunctions and eigenvalues of this linear prob¬ 
lem by tildes. In order to obtain a nonlinear branch 
of solutions (j)n which bifurcates from ^n{x) we require 
\imx^±oo 4^n{x) = 0. Thus all (j)n{x) can be considered 
real without loss of generality. Since Lq is Hermitian, all 
the eigenvalues jin are also real, what means the stabil¬ 
ity of the respective linear system (here we exclude the 
situation where zero is an eigenvalue of Lq). 

Turning now to e > 0 we observe that in the pres¬ 
ence of gain-and-loss distribution, stationary modes (if 
any) must have non-zero (hydrodynamic) current, i.e., x- 
dependent argument. Respectively we consider the con¬ 
struction of the T^T-symmetric potential and nonlinearity 
for the modes having the form 

(pn{x) = ^n{x,e)exp (^iej Vni^djj , (6) 

where the real function Vn{x) is the hydrodynamic veloc¬ 
ity and n stands here for the identification of the family 
bifurcating from the n-th linear eigenstate. Since this 
ansatz implies that the modulus of the nonlinear modes 
persists equal to the linear distribution, a possibility of 
constructing such modes is not obvious. In order to ad¬ 
dress this issue we substitute Eq. (j6]) in Eq. (j4]) and obtain 
relations linking the phase: 

^ (vn{x)\^n{x,€)\^'^ = 2W {x)\^n{x, e)f, (7) 

and the amplitude: 

^vl{x) +G{x)\(j)n{x,e)\‘^ + Vi{x) + e^V 2 {x) = Vn, (8) 

where we introduce the shift of the eigenvalue Vn deter¬ 
mined by /i — /in = 

The last equation has a free parameter e which leaves 
much freedom in constructing a solution. Although we do 
not require e to be small, we nevertheless restrict further 
consideration to two basic cases as follows: 

Case 1. The amplitude of (j)n{x) is independent on e, 
i.e. ^n{x,e) = ^n{x) and 

pi{x) + G{x)^l{x) + Vi{x) = Pn, V 2 {x) = 0 . ( 9 ) 

Case 2. The amplitude of (j)n{x) is proportional to e, 
i.e. (j)n{x,e) = e(j)n{x) and 

^vl{x) + Viix) = Pn, 


G{x)^‘l{x) + V2{x) = 0. (10) 
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The obtained system ([7]) with Eq. ([8j) [cf. Eq. (j9]) or 
Eq. (pT|) ] still contains freedom in the definition of the 
potentials and the hydrodynamic velocity. Therefore we 
impose further constraints on the part linear Vj{x) for 
j ^ 1, W{x) and nonlinear G{x) potentials, as well as 
the hydrodynamic velocity Vn{x)^ requiring them to be 
localized, i.e., Vj{x) ^ 0 for j > 1, ^ 0, G{x) 0, 

and Vn{x) ^ 0 at \x\ oo. Considering now Eq. ([8]) in 
the limit \x\ oc we readily obtain that = 0. Thus, 
in this case a nonlinear mode has not only the same form 
but also the same eigenvalue as its linear counterpart at 
e = 0 does. 

Thus we have obtained the profiles of the nonlinear 
modes in an exact analytical form t) = (l)n{x)e~^^^^ 
with 0n(^) being given by Eq. ([6j). Our next main task 
is to study the stability of the nonlinear modes. We do 
this numerically by two standard approaches. Eirst, we 
address the linear stability of a nonlinear mode 'ipnix, t) = 
(j)n{x)e~^^^^^ employing the ansatz 


= l^(j)nix) + g f{x)e + 5 *(a;)e*^** |e 

( 11 ) 

where p ^ 1, and f{x) and g{x) are the eigenfunctions 
of the linear eigenvalue problem: 


-K )\si) \9j ' ’ 

with 

Le = + Veix) + 2G^{x)\(pn{x)\'^ - fin- (13) 

The solution is linearly unstable if 6 has a non-zero imag¬ 
inary part, otherwise it is stable. 

Second, we test the stability by direct propagation us¬ 
ing an exact solution = (j)n{x)e~'^^^'^ with (j)n{x) 

given by Eq. (j6j) with a noise perturbation of order 1% 
of the initial amplitude \i^n{x^ 0)|, as the initial condition 
for Eq. ([1]). 


where Hn{x) = {d'^e ^^)/{dx^) is the Hermite 

polynomial with n = 0,l,2,..., and the eigenvalue is 

/i = /2n=n-h^, (16) 

in which we have = 0. 

We will also concentrate on the specific Gaussian form 
of the nonlinearity 

G{x) = 2ae-^^" (17) 

with the characteristic width Ij^/a with a > 0 and 
constant a (particularly, the nonlinearity is a constant 
for a = 0). This is a natural choice, for example, for 
the Bose-Einstein applications, where nonlinearity can be 
controlled through the optical Eeshbach resonance, with 
G(x) describing the profile of the laser beam (see, e.g.. 

Thus, the solution of the problem is reduced to two 
steps as follows. Eirst, given the gain-and-loss distribu¬ 
tion W{x) from Eq. ([7]) one obtains the hydrodynamic 
velocity: 

Mx)= r W{ 0 \MC,e)\^dC. (18) 

\(pn\X^ 6J I J — oo 

Second, the “correction” to the conservative part of 
the Hamiltonian is computed from Eqs. (P-([TQl). with 
G(x), Vn{x), and 0n(^) given by Eqs. (p!7|) . and ([IS]). 

In the meantime, formula m imposes the addi¬ 
tional constraint on the choice of the imaginary potential 
W (x) , which must ensure the existence of the integral in 
Eq. (fT8|) . This problem can be overcome if one again uses 
the inverse engineering, i.e. considers the hydrodynamic 
velocity given and finds W{x) from Eq. ([7|). 


IV. A SINGLE-WELL POTENTIAL 
A. PT-symmetry phases of the linear problem 


III. PERTURBED PT-SYMMETRIC LINEAR 
PARABOLIC POTENTIAL 

Below we concentrate on the physically relevant case 
of the parabolic (harmonic) potential Vo{x) = uj‘^x‘^j2. 
Without loss of generality, one can scale out the fre¬ 
quency uj making it one to yield: 

Vo{x) = (14) 

The families of the nonlinear modes of this potential at 
W{x) = X, Vj{x) = 0 (j = 1, 2) and constant nonlinearity 
G{x) =const were considered in Q. 

Now the profile of linear modes related to Eq. m is 
described by the Gauss-Hermite functions 

^n{x) = Hn{x)e~^^/‘^, (15) 


The simplest single-well potential is obtained by set¬ 
ting Vj{x) = 0 {j = 1,2), which corresponds to the 
Case 1 in Eq. [in the Case 2 in Eq. (nni) one has 
G{x)^‘^{x) = 0, which is the trivial solution]. Recalling 
that fin is given by Eq. (fT^ (i.e., Un = 0) and considering 
4>n{x) given by Eq. m we readily conclude that such a 
choice is possible only for attractive (focusing) nonlinear- 
ities cr < 0 (without loss of generality we set a = —1) for 
which the hydrodynamic velocity reads 

Vnix) = 

Now from Eq. (Q we can find the gain-and-loss distribu- 
tions (i.e. imaginary part of the potential) 

Wn{x) = [6nif„_i(x)-(a+3)a;if„(x)]e-('^i)"'/2_ ( 20 ) 

In order Wn{x) to be an odd function, i.e. to support 
the T^T-symmetry, we require n to be an even number: 
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FIG. 1. (Color online). Solid (red) and dashed (blue) curves 
indicate the lines of phase transitions for the T^T-symmetric 
potentials ([2^ and ([23l) with n = 0, 2. The unbroken (bro¬ 
ken) PT-symmetric phase is in the domain below (above) the 
phase breaking lines. Square-shaped (red) and circle-shaped 
(blue) curves indicate the lines of linear stability of solitons 
(IM)) with n = 0, 2. The stable (unstable) soliton is in the 
domain below (above) the linear stable lines. 

n = 0,2,4,...,. The functions Wn{x) with odd n are even 
and do not satisfy the condition of the T^T-symmetry. We 
do not consider them here but observe that for the odd n 
Eq. dU still obeys exactly localized solutions of the form 
m with Vo{x), ^n{x), G{x), Vn{x), and Wn{x) given by 
Eqs. ([15]) . ([T7|) . dH]), and ^0\i , respectively. 

We also observe that all the members of the Wn{x) 
family of potentials (i.e. the potentials corresponding 
to different n) are ^rt’o-parametric, i.e. are determined 
by the amplitude [it is given by e when substituted in 
Ue{x)] and by the internal parameter a. 

In order to establish domains of unbroken VT- 
symmetry phase of the linear T^T-symmetric potential 
Ue for different n we address the spectral problem 

Ln^{x) = Xn^{x), Ln = Lo-\-ieWn{x), (21) 

where Lq is given by Eq. An and 4^(x) are eigenval¬ 
ues and eigenfunctions, respectively. Since the discrete 
spectrum of a PT-symmetric potential is either real or 
appears in complex conjugated pairs, we conclude there 
exists a nonzero domain of the parameter e for which the 
PT-symmetry remains unbroken. 

The simplest potentials with nonzero complex parts 
are given by n = 0: 

U,{x) = — - ie{a + (22) 

and n = 2: 

[/,{a;) = y-2iea;[2(a + 3)a;2-(a + 15)]e-(“+i)^ /2. (23) 

In Fig. [T] we show the domains of broken and unbroken 
phase on the (a, e)-plane. The both curves in the figure 
grow with (a, what can be understood from the fact that 
growth of a corresponds to the shrinking of the gain-and- 
loss domains. 



FIG. 2. (color online). Real (a), (c) and imaginary (b), (d) 
parts of the eigenvalues An [see Eq. (EU] as functions of a for 
the potential Ea) atn = 0,e = 0.6 (the upper row) and for 
the potential ([23l) atn = 2, e = 0.18 (the lower row). 

In both cases, illustrated in Fig. |2l the spontaneous 
symmetry breaking occurs due to collision of the two 
lowest states as a decreases (what corresponds to the 
increase of the width of the gain-and-loss domains). All 
upper eigenstates (we checked numerically for the 6 low¬ 
est states) remain real. This can be understood from the 
fact that for the levels with large n the imaginary part 
of the potential represents a weak perturbation while the 
lowest levels are the most strongly deformed ones. We 
also notice that the instability is oscillatory (i.e. the two 
emergent complex eigenvalues have nonzero real part). 

B. Nonlinear modes in a single-well potential 

Now we turn to the nonlinear modes in the VT- 
symmetric potentials (j22j) and E3- The expression for 
the nonlinear modes is obtained form Eqs. m and m- 

(j)n{x) = if„(a;)e“*^/^exp[2zey . 

(24) 

Two examples of the modes are illustrated in Eig. [3] and 
their linear stability analysis is presented Eig.[TJ The fea¬ 
ture, most relevant to the present consideration, is dis¬ 
played by the domains between the solid line and squares 
for n = 0 and between the dashed line and circles for 
n = 2. In these domains the nonlinear modes [given 
by Eq. ([24]) with n = 0 and n = 2] are stable, while 
the respective linear PT-symmetric phases are broken, 
i.e. linear stability of the nonlinear modes is extended 
beyond the unbroken linear T^T-symmetric phase. Be- 
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FIG. 3. (Color online). One-hump nonlinear modes given by 
Eq. (fM)) with n = 0 and a = 1 for (a) e = 0.45 (real spectrum 
of the operator Lq , i.e. unbroken linear PT-symmetry) and 
(b) e = 0.7 (broken linear PT-symmetry). (c) Stable and 
(d) periodically varying propagation of the nonlinear modes 
described by Eq. (|24l) and corresponding to the weakly per¬ 
turbed initial conditions shown in (a) and (b), respectively. 


EIG. 4. (Color online). Two-hump nonlinear modes given by 
Eq. (|^ with n = 2 and a = 2 for (a) e = 0.02 (real spectrum 
of the operator L 2 , i.e. unbroken linear PT-symmetry) and 
(b) e = 0.12 (broken linear PT-symmetry). Both nonlinear 
modes are linearly stable [see Eig. [1] what are illustrated in 
panels (c) and (d) where the direct numerical simulations of 
Eq. ()24D are performed with weakly perturbed initial profiles 
(a) and (b). 


low we explore these domains in order to “draw” the 
mode along the branch bifurcating from the linearly sta¬ 
ble mode. Before that, however we show the check of 
stability by means of the direct propagation of the ini¬ 
tially stationary state Eq. (|24|) with a noise perturbation 
of order 1%. In Fig. [3] (c) we show stable propagation of 
the soliton for the parameters belonging to the domain of 
the unbroken linear T^T-symmetric phase of the operator 
I/O [defined in Eq. (j2T])], and to the linearly stable non¬ 
linear mode. In Fig. [3] (d) we illustrate the evolution of 
the mode where the linear T^T-symmetric phase is bro¬ 
ken and the nonlinear mode is linearly stable. In this last 
case we observe the oscillatory (breather-like behavior). 

Similarly, Figs. 0] (a) and (c) display the initial states 
and stable intensity evolution of a two-hump solitary 
wave [it is described by Eq. (|M|) with n = 2] for the 
parameters which guarantee both the real spectrum of 
the operator L 2 (unbroken T^T-symmetric phase) as well 
as the linear stability of the nonlinear mode. Meantime, 
Figs, m (b) and (d) show a two-hump soliton for the pa¬ 
rameters corresponding to broken liner T^T-symmetry 
but still stable the nonlinear mode. In both numeri¬ 
cal simulations we observed robustness of the nonlinear 
modes with respect to weak initial noise. 

Now we turn to the excitation of nonlinear modes by 
means of slow change of the control parameter e(t) which 
is now considered as a function of time. More specif¬ 
ically we consider simultaneous adiabatic switching-on 
the gain-and-loss distribution and the nonlinearity, mod¬ 


eled by [cf. Eq. ([T])] 


( 25 ) 

where the single-well potential V{x)^ nonlinearity G{x) 
and gain-and-loss distribution W(x) are given by 
Eqs. (p!4|) . (p!7|) and (|2Q|) . and 


e(t) 


0.2sin( )+0.45, 0 < t < 600 

V12OO; - (26) 

0.65, 600 < i < 1200. 


This choice of the final value of e is justified by the fact 
that the whole “trajectory” e(t), it is shown in Fig.[5][a), 
belongs to the parameter domain where the nonlinear 
mode is stable. 

Fig. [Sfb) exhibits the evolution of the solution '0(x,t) 
governed by Eqs. ([25]) and (j26|) subject to the initial con¬ 
dition given by (|M|) with n = 0, i.e., for the single-well 
potential Vo(x), nonlinearity G(x) and gain-and-loss dis¬ 
tribution W{){x) given by Eqs. ([T4|) . (fTTj) . and (j20|), re¬ 
spectively. One observes remarkably stable propagation 
with the increasing amplitude of nonlinear modes, which 
is driven from a nonlinear mode at the system parame¬ 
ters of unbroken linear T^T-symmetric phase, to the sta¬ 
ble nonlinear mode at the parameters where the liner 
PT-symmetric phase is broken (i.e., while e(t) is grow¬ 
ing there occurs the linear T^T-symmetry phase transi¬ 
tion and linear modes becomes unstable). 

Next we consider the excitation of the mode (|24|) with 
n = 2 described by Eq. ([25]) with the control parameter 
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FIG. 5. (Color online), (a) Time dependence of the pa¬ 
rameter e(t) (solid line) given by Eq. (|26l) . The squares and 
solid line (they correspond to the squares and dashed line 
in Fig. [T]) indicate the phase T^T-symmetry breaking phase 
transition in the linear case and loss of the stability of the 
nonlinear modes, respectively, (b) Intensity evolution of a 
nonlinear mode governed by Eq. (|25l) with the initial condi¬ 
tion 'ipo{x, t = 0) = 00(ic) given by Eq. (l2^ at e = 0.45 for the 
single-well potential Vo{x) ([T^ and gain-and-loss distribution 
Wo{x) (|20)) . Other parameters are n = 0 and a = 1. 



FIG. 6. (Color online) (a) Time dependence of the parameter 
e(t) (solid line) given by Eq. {|27|). The circles and dashed line 
are explained in Fig.[T] (b) Intensity evolution of a two-hump 
nonlinear mode of Eq. (|25ll with n = 2, and a = 2 subject to 
the initial condition 02(a:, t = 0) = 02(a;) given by Eq. (1^ 
at e = 0.02 for the single-well potential Vq{x) da and gain- 
and-loss distribution W 2 {x) (1^ . The dashed line shows the 
finally established amplitude of the mode which is different. 


adiabatically changing according to the law [see Fig. [6] 

(a)] 

{ 0.1 sin f +0.02, 0<t< 600, 

V 12 OO; - (27) 

0.12, 600 < t < 1200. 

The single-well potential Vb(x), nonlinearity G{x) and 
gain-and-loss distribution W 2 {x) given by Eqs. ([H]), (pTj) 
and (EOl), respectively. 

Like in the previous case we observe stable evolution of 
the nonlinear mode “prepared” in a system with unbro¬ 
ken T^T-symmetric phase and drawn to the system with 
broken T^T-symmetric phase. The amplitude of the non¬ 
linear mode grows with e(t) and remains unchanged after 
the control parameter reaches its final value [see Fig. [6] 

(b)]. 


V. A DOUBLE-WELL POTENTIAL 


A. T^T-symmetry phases of the linear problem 


Now we consider multi-well potential with 
Vi{x)V 2 {x) ^ 0, which corresponds to the Case 2 
in Eq. m- Recalling that fin is given by Eq. m and 
considering 0n(^) given by Eq. m we conclude that 
now the hydrodynamic velocity takes the form 

Vn{x) = (28) 

(i.e., it is now fixed having no free parameters). 

From Eqs. m and m it follows that the gain-and- 
loss distributions (i.e. imaginary potentials which are 
odd functions) numbered by an even number n, i.e. en¬ 
suring T^T-symmetry of the potential, and generated by 
the velocity field (j28j) are given by 


^n(x) = 


3nHn-i{x) - -xHn{x) 




(29) 


The real parts of the respective linear potential read 

V,nix) =-^H^^{x)e--\ (30) 

V 2 nix) =-2aHl{x)e-^‘^+^>\ (31) 

All the members of the Wn{x) family of potentials are 
two-parametric; they are determined by the amplitude 
given by e in U^{x)^ and by the internal parameter a 
introduced in Eq. m- 

Unlike in the Case 1 considered in the previous Section, 
for the case of a multi-well potential the parameter a [see 
Eqs. (p!7|) and (|3T]) ] need not to be non-positive and can 
acquire positive values, i.e. the nonlinearity can be either 
attractive or repulsive. Since however in this section we 
are interested in a double-well potentials, we concentrate 
on the “deformation” of the linear potential 

Ue{x) = ^ + e^Vio{x) + e‘^V 2 oix) + ieWo{x), (32) 
where 

Vio(a:) =^20(2:) =(33^) 

Woix) = -^xe-^"/\ (33b) 

[this corresponds to the chemical potential jx = flo = 1/2] 
and require a < — l/(4e^), the latter constraint ensuring 
the required double-well shape. 

Notice that at a = 0 the nonlinearity is a constant 
[see Eq. dni] while the real part of the potential U^{x) 
preserves a double-well shape. Nonlinear modes in this 
case were recently reported in (also, nonlinear 

modes in a slightly different double well potential with 
the imaginary part (j33b|) was investigated in [29| ). When 
a increases of the amplitude of the potential barrier be¬ 
tween the humps decreases. 
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FIG. 7. (Color online), (a) The line of the T^T-symmetry phase transition for the PT-symmetric potentials ()32l) . The unbroken 
(broken) T^T-symmetric phase is in the domain above (below) the phase breaking lines. Triangle-shaped curve indicates the 
border of linear stability domain of the mode (1351) with n = 0. Stable (unstable) solitons are in the domains above (below) the 
respective border lines. Real (b) and imaginary (c) parts of the eigenvalue A [see Eq. (IMl) ] as functions of a for the potential 
(l3^ at cr = —0.73. In all panels e = 0.8. 


As above we start by defining the domains of the un¬ 
broken and broken T^T-symmetric phase of the underline 
linear problem with n = 0. Respectively, we consider the 
linear spectral problem [cf. Eq. (|2T])] 

= X^, L = Lo + e^Vwix) + e^V 2 oix)+ieWoix), (34) 

where Lq is given by Eq. m, A and ^ = ^(x) are 
the eigenvalue and eigenfunction, respectively. For the 
present case it will be important that the nonlinearity 
can also be varied in time and its value affects the sta¬ 
bility of the mode. Indeed, in Fig. [3 (a) we show the 
domains of broken and unbroken phase at e = 0.8 on the 
(a, cr)-plane. Parametric dependence of the lowest eigen¬ 
values is shown in Figs. [3(b) and (c). Like in the case of 
one-hump potential we observe that the unbroken phase 
corresponds to relatively large a and the spontaneous 
symmetry breaking occurs as a decreases (what corre¬ 
sponds to the increase of the width of the gain-and-loss 
domains) due to the collision of the two lowest eigenval¬ 
ues. However now the broken phase corresponds to a 
limited interval of a and we observe the re-entered un¬ 
broken phase as a approaches zero. We also notice that 
the instability is oscillatory: the emergent complex eigen¬ 
values have nonzero real part. 


B. Nonlinear modes in a double-well potential 

Now we turn to the nonlinear modes in the VT- 
symmetric double-well potential ([32]), whose explicit ex¬ 
pression is obtained from Eqs. and (j28j): 

(pn{x) = J -(35) 

The results of the linear stability analysis of solution 
(|35|) are shown in Fig. [3^a). The feature most relevant 
for our consideration consist in the domain, now with 
respect to cr, where the linear T^T-symmetry is broken. 


while the nonlinear mode remains stable (in analogy with 
the case of the one-well potential). 

The stability of the nonlinear mode ([35]) with n = 0 
is also confirmed by the direct propagation with the per¬ 
turbation of initial profile, as illustrated in Fig. [51 In 
Fig. [8] (b) we show stable propagation of the soliton for 
the parameters belonging to the domain of the unbroken 
linear T^T-symmetric phase of the operator L [defined in 
Eq. (|34|)], and to the linearly stable nonlinear mode. In 
Fig. [8] (c) we illustrate the evolution of the mode where 
the linear PT-symmetric phase is broken, however the 
nonlinear mode is linearly stable. 

Now we turn to the excitation of nonlinear modes in 
the double-well potential. Since the potential V 2 o{x) and 
nonlinearity G{x) both contain <j, this parameter a can 
be exploited for managing (i.e. for excitation, in our case) 
of the nonlinear modes. To this end we consider a to be 
a function of t and address the adiabatic switch-on of the 
potential and the nonlinearity, governed by the model m 
which now is rewritten in the form 

idt'ip = + [Vo{x) + e^Vio{x) + e^V 2 o{x-, a{t)) 

+ieWo{x)]'tp + e^G'(x; (36) 


Here the double-well potential Vio{x) and gain-and-loss 
distribution Wo{x) are given by Eqs. (|33ap and (|33bp . 
respectively, and V 2 o{x; a{t)) G{x;a{t)) stand for V 2 o{x) 
and G{x) given by Eqs. (I33ap and ([TTP with a replaced 
by a{t). 

According the general idea described above, now we 
choose the adiabatic change of cr(t) in such a way that it 
assures that the system evolves from the domain of the 
unbroken T^T-symmetry of the underlying linear model 
to a broken phase, where however the nonlinear mode is 
linearly stable. More specifically, in the numerical simu¬ 
lations we exploit the dependence 


(j{t) 




—0.6 sin 



-1.25, 


0.65, 0 < t < 600, 

(37) 

600 < i < 1200. 
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FIG. 8. (Color online). The same nonlinear modes given by Eq. (l35l) with n = 0, a = 4, e = 0.8 for both a = —0.65 
(real spectrum of the operator L, i.e. unbroken linear T^T-symmetry) and a — —1.25 (broken linear PT-symmetry) (a). 
Both nonlinear modes are linearly stable [see Fig. [Zl[a)] what are illustrated in panels (b) and (c) where the direct numerical 
simulations of Eq. (EH) are performed with weakly perturbed initial profile (a) for different parameters a = —0.65 and a = —1.25, 
respectively. 



FIG. 9. (Color online), (a) Time dependence of the param¬ 
eter (j(t) given by Eq. (|37l) (solid line). Upper dashed line 
and triangles indicate boundaries the stability of the nonlin¬ 
ear mode and the unbroken PT-symmetric phase (in both 
cases above the respective curve), (b) Intensity evolution of 
nonlinear mode of Eq. (I36p with the initial condition given by 
Eq. (l35)) for the double-well potential with n = 0, a = 4, and 
e = 0.8. 


The dependence cr(t) is illustrated in Fig. [9] (a). 

Fig. [9] (b) exhibits the excitation of a nonlinear mode 
described by Eq. (|36|) for the double-well potential, i.e. 
Vo{x), Vio{x), V 2 o{x; (j{t)), nonlinearity G{x]a{t)) and 
gain-and-loss distribution Wq (x) determined by Eq. (l37|) 
for the varying parameter The initial condition 

in these simulations is taken the form dm with n = 0. 
In the figure we again observe the stable evolution of 
the nonlinear mode between the initial and final shapes 
of the potential (notice that the change of the mode 
amplitude is relatively small and not clearly visible on 
the scale of the figure). 


VI. DISCUSSION AND CONCLUSIONS 

In the present paper we reported several new 
branches of nonlinear modes described by the nonlinear 
Schrodinger equation with the PT-symmetric single-well 
and double-well potentials and Hermite-Gaussian distri¬ 


butions of the gains and losses. The reported solutions 
are two-parametric, each of the parameters defining a 
branch of the solutions. All the considered branches bi¬ 
furcate from the modes of the respective linear poten¬ 
tials. A peculiarity of the reported modes consists in 
their stability properties: their stability in the param¬ 
eter space extends beyond the domains of the stability 
of the respective limits, i.e. beyond the domain of the 
unbroken PT-symmetric phase of the underlying linear 
problem. This suggests a possibility of how a stable non¬ 
linear mode can be excited in a T^T-symmetric system 
with broken T^T-symmetry. The method is based on the 
drawing the mode by adiabatic change of one of the con¬ 
trol parameters along the branch from the domain where 
linear stability is verified and the PT-symmetric phase 
is unbroken, i.e. where the mode can be excited by one 
of the conventional methods, to the domain where the 
PT-symmetry is broken but the nonlinear stability still 
persists. 

While the described approach is straightforward math¬ 
ematically it leaves an open problem on managing gains 
and losses by a single parameter, while in many cases the 
physical phenomena responsible for gains and losses have 
different nature. Speaking more generally, to implement 
the suggested scheme in practice one needn’t to hold an 
integrable model but one does need to preserve balance 
between the gains and losses. This last goal be achieved 
in at least one of the ways as follows. First, consider¬ 
ing PT-symmetric profiles cerated in mixtures of atomic 
gases [I 3 5 oiie can modify both active and lossy domains 
by a single parameter, which is the intensity or wave¬ 
length of the control field, or alternatively by varying in 
space the properties of cladding of the atomic cell. The 
latter affects the whole spectrum of the underlying lin¬ 
ear problem, i.e. both real and imaginary parts of the 
refractive index. On the other hand when considering 
PT-symmetric double-well potentials in the BEG prob¬ 
lems, to ensure the balance between varying gains and 
losses one has to perform simultaneous loading atoms to 
one of the wells (either using atomic laser or tunnel¬ 
ing from a neighbor potential well, as suggested in [T^) 
and eliminating atoms (using, for example, ionization by 
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external beam [33| or leakage of atoms through tunnel¬ 
ing [ill). 

Finally, returning to the considered exact model, the 
method was tested using the found exact solutions in 
the form of one- and two-hump modes supported by the 
one-well T^T-symmetric potential and by the repulsive 
nonlinearity, as well as for the modes in a double-well 
potential and attractive nonlinearity. In all the cases we 
observed the stable evolution of nonlinear modes, thus 
supporting practical feasibility of the approach. 
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